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We review the salient features of _B-meson physics, with particular emphasis 
on the measurements carried out at the B-factories and Tevatron, theoreti- 
cal progress in understanding these measurements in the context of the stan- 
dard model, and anticipation at the LHC. Topics discussed specifically are the 
current status of the Cabibbo-Kobayashi-Maskawa matrix, the CP-violating 
phases, rare radiative and semileptonic decays, and some selected non-leptonic 
two-body decays of the B mesons. 

1. Introduction 

The two B-meson factories operating at the KEK and SLAC e + e~ storage 
rings have outperformed their projected luminosities and have produced 
a wealth of data on the process e + e~ — * T(45) — * BB, with subsequent 
weak decays of the B and B mesons. The exclusivity of the final state, the 
asymmetric beam energies and the large statistics collected at the T(4S) 
(O(l0 9 )BB events) have led to a number of impressive and quantitative 
results, which include, among other measurements, precise determination 
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of the weak mixing (Cabibbo-Kobayashi-Maskawa CKM ' ) matrix ele- 
ments | Vcb | and \V u b\ ^, a large number of CP-violating asymmetries in 
exclusive decays of the B mesons ^, and rare radiative and semileptonic de- 
cays (implying flavor changing neutral current FCNC transitions) involving 
the B° and B + mesons and their charge conjugates ^. Recently, Belle col- 
laboration have published a number of interesting results on the decays of 
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the B s -meson highlighted by the first measurement of the radiative penguin 
decay B s — > 4>j ^ . This is the SU(3)-counterpart to the much celebrated 
decay B — > K*j, first measured by CLEO and further consolidated by the 
BABAR and BELLE collaborations 1 . 

The two B-meson factories are not the only players currently active in 
the study of £>-meson physics. In fact, 6-physics is the main beneficiary of 
the upgrade program at the high energy proton-antiproton collider Tevatron 
at Fermilab. Since this upgrade, physics news from Fermilab are dominated 
by the achievements in the &-quark sector. We have seen very impressive 
and seminal results in this field by the two Fermilab experiments CDF and 
DO, including the landmark measurement of the B® - S^-mixing induced 
mass difference AMb 3 ^' ^ and the first result on direct CP violation in 
the i^-meson sector Acp(B® — ► K + tt~) ^ These experiments have 
established (if any proof was needed) that cutting edge flavor physics is 
done also at hadron machines. In all likelihood, this success story of the 
hadron machines will be set forth at the Large Hadron Collider (LHC), 
commissioned to operate in 2008 at CERN. In particular, the experiment 
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LHCb , dedicated to precision measurements of the physics of the entire 
family of B-mesons (B°, B + , B®, their charged conjugated states, 
and excited states) and Ab baryons (including the entire baryon spectrum 
containing at least one 6-quark), will greatly broaden our knowledge of both 
the spectroscopic and dynamical aspects of 6-physics. The other two LHC 
experiments ATLAS and CMS will also contribute to b physics. 

In this talk, I will briefly review the highlights of the measurements al- 
ready accomplished, making contacts with theoretical expectations in the 
standard model (SM). The topics discussed are: (i) an update of the CKM 
matrix, with emphasis on the matrix elements in the third row and the 
third column of Vqkm from direct decays (which determine \V c b\, \ V u b\ and 
|Vffc|), and from induced transitions involving the mass differences AMs d , 
AMb s and the electromagnetic penguin decays b — > (s, d)-f (which deter- 
mine the matrix elements \Vtd\ and |Vt a |) (ii) measurements of the CP- 
violating phases a, (3 and 7 from exclusive i?-decays, and (iii) radiative, 
semileptonic and leptonic rare -B-decays, a field which has received a lot of 
theoretical attention both in the context of the SM and in extensions of it. 
Much more detailed discussions of these topics can be found in the original 
literature cited below, in the reviews, such as the Review of Particle Prop- 
erties by the Particle Data Group and in the proceedings of the topical 
workshops on flavor physics 
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2. Status of the CKM Matrix 

The CKM matrix is written below in the Wolfenstein parameterization ^ 
in terms of the four parameters A, A, p, r/: 



( 



V, 



CKM 



1-^ 

-A(l + iA 2 X 4 rj) 



X AX 3 (p~ir]^ 



\-\X z 



^AX 3 {1-p- vq) -AX 2 (l + iX 2 r^j 



(1) 



Anticipating precision data, a perturbatively improved version ^ of the 
Wolfenstein parameterization will be used below with p = p(l — A 2 /2), fj — 
V (l-X 2 /2). 




Fig. 1. The unitarity triangle with unit base in the p - r\ plane. 



Unitarity of the CKM matrix implies six relations, of which the one 
resulting from the equation V u dV* b + V c dV* b + V t dV t * b — is the principal 
focus of the current experiments in B-decays. This is a triangle relation in 
the complex plane (i.e. p~fj space), and the three angles of this triangle are 
called a, [3 and 7, with the BELLE convention being <f>i = /?, (f>2 = a and 
4>3 = 7 ( see Fig.l). The unitarity relation in discussion can also be written 
as 
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where Rb = 



Vab 

and the 



= V ' P 2 + r f an d Rt = x 
p) 2 + f] 2 . Thus, precise determination of |V c b|, \V u b\ and \Vtd 
three CP-violating phases a, (3, 7 is crucial in testing the CKM paradigm. 

In this section, we mainly review the current status of the CKM matrix 
elements in the third row and the third column, namely V u b, V c b, Vtd, Vts 
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and Vtb- The other four CKM matrix elements, V u d, V us , V c d and V cs are 
discussed in ^ ^, from where we quote the current best measurements. 

• V u d is obtained from the nuclear /3-decay of the (super-allowed) 
0+ -» (9+ transitions 17 < 18 , yielding: = 0.97377(27). 

• V^ s is obtained from the JQ 3 decays and chiral perturbation theory 
yielding V us = 0.22535 ± 0.00116. In getting this, the Leutwyler- 
Roos value for the form factor /+(0) = 0.961 ± 0.008 ^ has been 
used. This is to be contrasted with the value, V us = A = 0.2265(7), 
obtained from the unitarity condition on the first row of the CKM 
matrix and the input value for V u d given above. They arc in good 
agreement with each other, providing a precision test of the uni- 
tarity of the first row of Vckm- 

• V cs is obtained from the decay D — ► Kev e and the Lattice-QCD 
result for the D — > K form factor yielding \V CS \ — 0.996 ± 
0.008 ± 0.015 ± 0.104, where the first two errors are experimental 
and the third (dominant) error is from theory. 

• Despite impressive progress in the measurement of the Cabibbo- 
suppressed D-meson decays and Lattice-QCD, the determination 
of V c d is still dominated by the neutrino-nucleon production of the 
charm quark, yielding: \V c d\ = 0.230 ± 0.011 Within errors, 
\V c d\ — \ V US \, as anticipated from eq. (1). 



2.1. Current determinations of \V c b\ and \ V u b\ 

Determinations of \V c b\ are based on the semilcptonic decay b — > clvt.. 
This transition can be measured either inclusively through the process 
B — > X c £i/£, where X c is a hadronic state with a net c-quantum num- 
ber, or exclusively, such as the decays B — > (D,D*)£i/£. In either case, 
intimate knowledge of QCD is required to go from the partonic process to 
the hadronic states. The fact that irib 3> Aqcd has led to novel applications 
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of QCD in which heavy quark expansion (HQE) plays a central role and 
the underlying theory is termed as HQET. 

2.1.1. \V c b\ from the decays B — > X c I ut 

Concentrating first on the inclusive decays, the semileptonic decay rate can 
be calculated as a power series 

r = r + — r! + ^2r 2 + ^3r 3 + ..., (3) 

m b ml ml 
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where each I\ is a perturbation series in a s (rrib), the QCD coupling constant 
at the scale mj. Here T is the decay width of a free 6-quark, which gives 
the parton-model result. The coefficient of the leading power correction Ti 
is absent and the effect of the 1/m^ correction is collected in T 2 , which 
can be expressed in terms of two non-perturbative parameters called Ai - 
kinetic energy of the 6-quark - and A2 - its chromomagnetic moment. These 
quantities, also called fjl and respectively, in the literature, are defined 

9^ 9/1 9^ 0f\ 

in terms of the following matrix elements ■ ■ ■ °: 



where is the covariant derivative and heavy quark fields are character- 
ized by the 4- velocity, v. At C(Aq CD /m^), six new matrix elements enter 
in T 3 , usually denoted by pi j2 and 71,2,3,4- 

Data have been analyzed in the theoretical accuracy in which correc- 
tions up to O(a^Po), 0(a s kQcv/mb) and 0(Aq CD /t7?|) are taken into 
account z '' , with /3 being the lowest order coefficient in the expansion 

of the QCD /3-function. In addition to the other parameters, a quark mass 

97 

scheme has to be specified. Bauer et al. ' have carried out a comprehen- 
sive study of the scheme dependence using five quark mass schemes: 15, 
PS, MS, kinematic, and the pole mass. 

To extract the value of \V c b\ and other fit parameters, three differ- 
ent distributions, namely the charged lepton energy spectrum and the 
hadronic invariant mass spectrum in B — > X c £Dg, and the photon en- 
ergy spectrum in B — > X s -f have been studied. Theoretical analyses are 
carried out in terms of the moments. Defining the integral R n (E cut , fi) = 
J E dEi(Ei — fj,) n dT /dE(, where E cut is a lower cut on the charged lep- 
ton energy, moments of the lepton energy spectrum are given by (Ef) = 
Rn(E cu t,0)/Ro{E C ut,0)- For the B — > X c lvt hadronic invariant mass spec- 
trum, the moments are defined likewise with the cutoff E cut . Analyses of 
the data along these lines have been presented by a large number of experi- 
ments. For a summary and references to the original literature, see HFAG ^ . 

The global fit of the world data in the so-called IS'-scheme for the b- 
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quark mass undertaken by Bauer et al. leads to the following fit values 
for \V c b\ and m\ s : 



2M B X 1 = (B(v)\Q v (iD) 2 Q v \B(v)), 
6M B X 2 = (B(v)\Q v a^{iD^iD»}Q v \B(v)) , 



(4) 



\V cb \ = (41.78 ± 0.30 ± 0.08 TB ) x 10 
ml s = (4.701 ± 0.030) GeV. 
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(5) 



The BABAR collaboration have studied the dependence of the lepton and 
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hadron moments on the cutoff -E C ut and compared their measurements with 
the theoretical calculation by Gambino and Uraltsev ^° using the so-called 
kinematic scheme for the 6-quark mass m\™{p), renormalized at the scale 
/i = 1 GeV. Agreement between experiment and theory allows to determine 
the fit parameters in this scheme with the results . 

\V cb \ = (41.91 ± 0.l9 exp ± 0.28 H qe ± 0.59r al ) x 10~ 3 , 

m b (lGeV) = (4.613 ± 0.022 exp ± 0.027/^) GeV, 

m c (lGeV) = (1.187 ± Q.033 exp ± O.OAOhqe) GeV . (6) 

The two analyses (5) and (6) are in excellent agreement with each other. 
The achieved accuracy 8\ V c b\/\ V c b\ — 2% is impressive, and the precision 
on m& is also remarkable, Snib/mb = O(10~ 3 ), with a similar precision 
obtained on the mass difference ra b — m c . 

2.2. \V cb \ from B — > (D,D*)£u e decays 

The classic application of HQET in heavy — ► heavy decays is in the decay 
B — » D*£i/£. The differential distribution in the variable oj(= Vb-Vd»), 
where vb{vd') is the four-velocity of the £?(D*)-meson, is given by 

£ = l^l 2 - ™z>) 2 - 1) 1/2 |^)| 2 , 

where ^(w) is a phase space factor with ( 1 ) = 1, and T(lo) is the Isgur- 
Wise (IW) function ^9 with the normalization at the symmetry point 
!F(1) — 1. Leading Aqcd/w;, corrections in .F(l) are absent due to Luke's 
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theorem . Theoretical issues are the precise determination of the second 
order power correction to T(u = 1), the slope p 2 and the curvature c of the 
IW-function: 

T{u>) = ^(1) [1 - p 2 (to - 1) + c (w - l) 2 + ...] . 

Bounds on p 2 have been obtained by Bjorken ^1 an d Uraltsev which 
can be combined to yield p 2 > 3/4. Likewise, bounds on the second (and 
higher) derivatives of the IW-function have been worked out by the Orsay 
group 33 ; yielding c > 15/32. Apparently, the data sets used in this analysis 
differ significantly from experiment to experiment, resulting in considerable 
dispersion in the values of jF(l)|T4b| and p 2 and hence in a large x 2 of the 
combined fit, summarized by HFAG : 

J"(l)|V cfc | = (35.89 ± 0.56) x 1CT 3 , (7) 
p 2 = 1.23 ± 0.05 (x 2 = 37.8/17; CL = 0.026). 
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To convert this into a value of \V c b\, we need to know T(X). In terms of the 
perturbative (QED and QCD) and non-perturbative (leading 5i/ m 2 and 
sub-leading <$i/ m 3) corrections, can be expressed as follows: 

^(1) = VA [l + £l/m 2 +h/ m *\ , (8) 

where t\a is the perturbative renormalization of the IW-function, known 
in the meanwhile to three loops . One- and two-loop corrections yield 
t\a — 0.933 and the 0{a 3 s ) contribution amounts to ij^ = —0.005. A 
Lattice-QCD calculation in the quenched approximation yields ^ = 
0.9191q;q35, which is now being reevaluated with dynamical quarks. Taking 
into account this theoretical input, the value quoted at the Lepton-Photon- 
2007 Symposium is 3 : 

(39.1 ± 0.65 exp ± 1.4 theo ) x 10- 3 . (9) 

The resulting value of \V c b\ is in excellent agreement with the ones given in 
(5) and (6) obtained from the inclusive decays. 

2.3. \V u b\ from the decays B — > X u £f£ 

HQET techniques allow to calculate the inclusive decay rate B — ► X u £vi 
rather accurately. However, the experimental problem in measuring this 
transition lies in the huge background from the dominant decays B — ► 
XJLvt which can be brought under control only through severe cuts on 
the kinematics. For example, these cuts are imposed on the lepton energy, 
demanding Eg > {m B — m 2 D ) /2ms , and/or the momentum transfer to the 
lepton pair q 2 restricting it below a threshold value q 2 < q^^ and/or 
the hadron mass recoiling against the leptons, which is required to satisfy 
mx < mjj. With these cuts, the phase space of the decay B — > X u tvt is 
greatly reduced. A bigger problem is encountered in the end-point region 
(also called the shape function region), where the leading power correction 
is no longer 1/m 2 but rather l/m&AQCD, slowing the convergence of the 
expansion. Moreover, in the region of energetic leptons with low invariant 
mass hadronic states, the differential rate is sensitive to the details of the 
shape function f(k+) where k + = k + k 3 with <~ O(Aqcd)- 

The need to know f(k + ) can be circumvented to a large extent by doing 
a combined analysis of the data on B — > X u lvt and B — ► X s j. Using the 
operator product expansion (OPE) to calculate the photon energy spec- 
trum in the inclusive decay B — > X s ~{, the leading terms in the spectrum 
(neglecting the bremsstrahlung corrections) can be re-summed into a shape 
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function 



37. 



£_£^ivs.iai"i<ri"/<i-.>. no) 

where x = and C5 S is an effective Wilson coefficient, characterizing the 

nib ' ' 1 & 

strength of the electromagnetic dipole operator. In the leading order, Eg- 
and Mx„-spectra in B — ► X u £vt are also governed by /(a;). Thus, /(a;) can 
be measured in B — > X,7 and used in the analysis of data in £? — ► X u £i/g. 
Following this argument, a useful relation emerges 40 

Kb '-^Vfi'^Va+w), (id 



m B /2 



r u (E c )= / dE, 



2 f mB/2 ,„ ,„ „ , dT s 



where 



T S (E C ) = — / dE^ -E c )—^, (12) 
m b J Ec d£/ 7 

and S(E C ) incorporates the sub-leading terms in C(Aqcd/ti6), which can 
only be modeled at present. In addition, there are perturbative corrections 
to the spectra and in the relation (11) ^7, 41, 38^ 

The other strategy is to extend the measurements of the inclusive decay 
B — > X u l + V£ to the kinematic regions which are reachable by the decay 
B — > X c £ + i>e, thus enlarging the region where the light-cone momentum 
component satisfies P+ = — 3> Aqcd, obviating the need to know 
the shape function. Both these methods have been used in the determina- 
tion of | V u b |, summarized below. 

• Determination of \V u b\ using the combined cuts on the variables m x 
and q 2 following the suggestion by Bauer, Ligeti and Luke (BLL) ^ . Taking 
the HQE parameter input from the analysis of the decay B — > X c £i/£ and 
B -> X s7 , HFAG 7 quotes an average value |Kb|(BLL) = (4.83 ± 0.24 ± 
0.37) x 10~ 3 , using the b-quark mass m 6 (lS) = (4.70 ± 0.03) GeV in the 
so-called lS'-scheme. 

• Determination of \V u b\ from a fully differential decay rate for 
B — > X u £vi based on the soft collinear effective theory (SCET) tech- 
niques 43, 44> 46^ hereinafter called the Bosch-Lange-Neubert-Paz 
(BLNP) approach 47 The three independent kinematic variables arc cho- 
sen to be: P e = M B -2E e , P_ = E X + \P X \, and P+ = E X -\P X \, where P± 
are the light-cone components of the hadronic final-state momentum along 
the jet direction, E x is the jet energy, P x is the jet momentum, and Eg is 
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the charged-lepton energy. In terms of these variables, the triple differential 
distribution is: 



^kdF e = ^TPM^ - P + ) {(P- + P + )(M B -P- + Pe~ P + ) Fl 
+ (Mb - P-)(P- - P+)F 2 + (P- - Pe){Pt - P+)F 3 } . (13) 

The "structure functions" F t can be expressed as product of the hard (per- 
turbatively calculable) coefficient and a jet function, which are convoluted 
with the soft light-cone distribution functions, the shape functions of the 
B meson. SCET allows to separate the two scales here, namely fih ~ mb 
and Hi ~ \Jm\, Aqcd and enables to sum large logarithms involving the 
two scales jj,h and /Ltj. The dependence on the sublcading shape functions is 
studied by taking several models. Fixing the HQE parameters (in the so- 
called shape function scheme) to the values m&(SF) = (4.63 ± 0.06) GeV, 
/Li^.(SF) = (0.18±0.06) GeV 2 , and the exponential form for the shape func- 
tion, HFAG 7 quotes \V ub \ = (4.31 ± 0.17 ± 0.35) x 10~ 3 in the BLNP 
approach. 

• Determination of \V u b\ using the so-called Dressed Gluon Exponen- 
tiation (DGE) advocated by Andersen and Gardi . The basic assump- 
tion of this approach is that properly defined quark distribution in an on- 
shell heavy quark provides a good approximation to the distribution in 
the meson. The problematic small hadronic mass Mx u region in the decay 
B — > X u l + i>i is characterized by a large hierarchy in the ratio p\ jp~ <C 1 
involving the partonic light-cone coordinates pj < p 3 < mj. Defining the 
moments n with respect to the powers of 1 — Pj^/pJ , the region of small p^ 
is probed by high moments n — > oo, giving rise to the Sudakov logarithms 
Inn. Infrared sensitivity appears in the moment-space Sudakov exponents 
through infrared renormalons, leading to the divergence of the (higher or- 
der perturbative) series in the Sudakov exponent. However, the leading 
infrared renormalon ambiguity cancels exactly against the pole-mass renor- 
malon ambiguity, which enters as a kinematic factor in the differential 

decay width, making the on-shcll perturbative calculation directly appli- 
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cable for the phenomenology This approach, applied to the data yields 
| Kb | = (4.34±0.16±0.25)xl0- 3 for the input value rn^(fh b ) = (4.20±0.07) 
GeV. 

Thus, the three theoretical approaches used in the determination of \V u b\ 
from the inclusive decays B — > X u £u£ give very consistent values. 
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2.4. \V u b\ from exclusive decays 

\V u b\ has also been determined from the exclusive decays B — > 7r^^. The- 
oretical accuracy is limited by the imprecise knowledge of the form fac- 
tors. A number of theoretical techniques has been used to determine them. 
These include, among others, Light-cone QCD sum rules ^ , Quenched- and 
Unquenched-Lattice QCD simulations 51, 52 jjpAG quotes the follow- 
ing extracted values of |Kb| using the full g 2 -range for the form factors (in 
units of 1CT 3 ) 7 : 

|Kb| = 3.43±0.10lg;S (LCSR: Ball-Zwicky 49 ); 
| Kb | = 3.17 ±0.10^;2 (Lattice: HPQCD 50 ); 
|Kb| = 3.82 ± 0.12±g;|| (Lattice: FNAL 51 ); 
|Kb| = 3.61 ± O.lliJ;^ (Lattice: APE 52 ). 

These values of |K&| from the exclusive decays (typically 3.5 x 10~ 3 ) are 
smaller than the ones extracted from the inclusive measurements (typically 
4.3 x 1CP 3 ). This mismatch is a mild irritant, waiting for its resolution 
through progress in Lattice-QCD and more B-factory data. 

To reduce the form-factor related uncertainties in extracting |Kb| from 
exclusive decays B — ► {it,p)lvi, input from the rare _B-decays B — > 
(K, K*)£ + £~ and HQET may be helpful. A proposal along these lines is 
the so-called Grinstein's double ratio which would determine |Kb|/|14bKsl 
from the end-point region of exclusive rare i?-meson decays To carry out 
this program one has to measure four distributions in the decays B — > p£f£, 
B -> K*£+£-, and D -» (p, K*)lv t . With the help of this data and HQET, 
the ratio of the CKM factors |Kb|/|VtbKsl can be determined through the 
double ratio 

T(B -> gfr) r(D -> 
T{B ^ K*£+£-) T(D^plv) ' 1 j 

At the -B factories, one expects enough data on these decays to allow a 10% 
determination of |K&| from exclusive decays. 

3. Status of the Third Row of Vckm 

FCNC transitions b — > s and b — > d (as well as s — > d) give information 
on the third row of the CKM matrix Vckm and allow to search for physics 
beyond the SM. In the SM these transitions are generally dominated by 
the (virtual) top quark contributions giving rise to the dependence on the 
matrix elements |V^Vt s | (for b — > s transitions) and |V^Vt<i| (for b — > d 
transitions). Of these, only the matrix element |Vth| has been measured by 
a tree process t — > Wb at the Tevatron through the ratio 
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Bjt^wb) \v tb f 

nth ~ B(t - Wq) \V td f + \V ts \i + \V tb \ 2 ■ [ °> 

The current measurements yield R tb — 1.12to;ig(stat)l^i3( s 2/ s ^)) yield- 
ing R tb > 0.61 (at 95% C.L.), which in turn gives \V tb \ > 0.78 at 95% C.L.. 

Precision on the direct measurement of |Vtf>| is rather modest (unitarity 
gives \Vtb\ — 0.9992.), which will be greatly improved, in particular, at a 
Linear Collider but also at the LHC. The corresponding measurements 
of | Vt s | and \ V t d\ from the tree processes are not on the cards. These matrix 
elements are determined by (loop) induced processes discussed below. 

3.1. Status of \V td \ 

The current best measurement of \V t d\ comes from AM^, the mass differ- 
ence between the two mass eigenstates of the B^ - B^ complex. This has 
been measured in a number of experiments and is known to an accuracy of 
~ 1%; the current world average is 7 AM Bd = 0.507 ± 0.005 (ps) _1 . 

In the SM, AM Bd and its counterpart AM Bsl the mass difference in 
the .B° - B® system, are calculated by box diagrams, dominated by the 
Wt loop. Since (M^,m t ) 3> m/,, AM Bd is governed by the short-distance 
physics. The expression for AM B taking into account the perturbative- 
QCD corrections reads as follows ™ 

AM b, = ^ Vb \V td V t * b \ 2 M Bd {fl d B Bd ) Ml S (x t ) . (16) 

The quantity t) b is the next-to- leading log (NLL) perturbative QCD renor- 
malization of the matrix element of the (|A£?| = 2, AQ — 0) four-quark 
operator, whose value is fj B = 0.55 ± 0.01 x t = m^/M^ and 

So(xt) — Xtji{xt) is an Inami-Lim function with 

,..191 31 3 s 2 lni 

h{x) = 4 + 4 (I" V) 2{l^xY 2 (1 = Icf ■ 

The quantity f Bd B Bd enters through the hadronic matrix element of the 
four-quark box operator, defined as: 

<5°l(Ml - 75)<Z)X°> = \f B B Bq Ml q , (18) 

with B q = B d or B s . With AM Bd and i) B known to a high accuracy, and the 
current value of the top quark mass, defined in the MS scheme, fht(mt) = 
(162.3 ±2.2) GcV, leading to S (x t ) = 2.29(5), i.e., 6S (x t )/ S (x t ) ~ 2.5%, 
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the combined error from these sources is a few percent. This is completely 
negligible in comparison with the current theoretical uncertainty on the 



matrix element f Bd y B Bd . For example, 0(a s )-improved calculations in the 

QCD sum rule approach yield f Bd = (210±19) McV 59 and f Bd = (206±20J 
MeV ^O, whereas B Bd in the MS scheme in this approach is estimated as ^ 
B Bd = 1 to within 10%, yielding for the renormalization group invariant 



quantity B Bd ~ 1.46, and an accuracy of about ±15% on f Bd y B Bd . Lattice 



calculations for f Bd y B Bd are uncertain due to the chiral extrapolation. 
Taking this into account, the current unquenched lattice QCD calculations 
from the JLQCD and HPQCD Collaborations yield 62 - 63 



fB d \f- 



B Bd = (244 ±11 ±24) MeV, (19) 

where the first error is statistical, and the second combines the uncertainty 
from chiral extrapolation and other systematic errors. Using \Vtb\ = 1, this 
yields 64 

\Vtd\ = (7.4 ±0.8) x 1(T 3 . (20) 



3.2. Present status of \Vt s \ 

The quantity which currently provides the best determination of |Vt s | is 
AM Bg , measured by DO and CDF, with the CDF measurement being more 
precise AM Bs = (17.77 ± 0.10± 0.07) (ps)" 1 11 . The expression for AM Bs 
in the SM can be obtained from the one for AM Bd (16) by the replacements: 
M Bd -> M Bs , {P Bd B Bd ) -» (f Bs B Bs ), and, most importantly, V td -» V ts . 

Using the recent Latticc-QCD result f Bs \[B~ B ~ s = (281 ± 21) McV from 

the HPQCD collaboration ^ allows to determine |Vt«| from AM Bs with 
a precision of about 10%, comparable to the precision on \V t d\ in eq. (20). 
This is not competitive with the indirect estimate of this matrix element 
from the CKM unitarity, which yields V ts ~ -V cb = -4.1(1) x 10~ 2 . 

The compatibility of the SM with the measured value of AM Bs is usually 
tested by taking the value of Vts from the unitarity fits, the measured value 

of m t , and the Lattice-QCD value for f Bg \J B Bs . Following this reasoning, 
the Lattice- HPQCD collaboration estimates AM Bs {$M) = (20.3±3.0±0.8) 
(ps) _1 ^ . The corresponding estimate from the UTfit is AM Bs (SM) = 
(20.9 ± 2.6) (ps)- 1 66 and the CKMfitter yields AM Bg (SM) = (21.7±|;|) 
(ps) _1 ^ . Thus, the experimental measurement of AM Bs is typically about 
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la below the SM estimates, with AM s (expt)/AM s (SM) = 0.88 ± 0.13 65 , 
0.85 ± 0.10 66 and 0.88 ± 0.20 67 . 

Possible implications of the CDF measurement of AM Bs have been 
studied in several papers ^8, 69, 70, 71, 66 p erna p Sj it i s to the point to 
mention here that a value of AMb s — 18 (ps) _1 was also hinted by the 
LEP data with the central value being 17.7 (ps) _1 . Anticipating this, the 

consequences of an eventual measurement of AMr around this value were 

79 

worked out some eight years ago ' for the parameters of the SM and 
the minimal flavor- violating supersymmctry .It was emphasized that a 
measurement of AMr around this value would confirm the SM. The CDF 
measurement of AMb s is very precise, and our knowledge of the CKM 
parameters and the non-perturbative quantities has in the meanwhile also 
improved, as discussed in this report. However, the bottom line remains the 
same, namely that the SM has passed this crucial test comfortably. 

The ratio of the mass differences AM^/AM^, now measured very 
precisely, can be used to constrain the CKM ratio | V t d\/\ V ts \ using the SM 
relation 74 : 

AM Bs = Mg. \V t lVts\ 2 . 
AM Bd $ M Bd \V* b V td ? ' 1 ' 



where £ = Jb s y Bb s / fB d \/ Bs d - Theoretical uncertainty in £ in the Lattice 



QCD approach is arguably smaller compared to the one in fs s y Bb s , as in 
the SU(3) limit £ = 1, and the uncertainty is actually in the SU(3)-breaking 
corrections. Current estimate in the unquenched lattice calculations of £ 

is 64 f = i.2ii°-°35, which y idds 11 

\Vtd/V ta \ = 0.2060 ± 0.0007(cxp)t0;^(th). (22) 

This is by far the best measurement of this CKM ratio, and it provides a 
non-trivial constraint on the allowed profile of the unitarity triangle. Com- 
bining eqs. (22) and (20) yields |V ts | = (36 ± 4) x 10~ 3 with the error 
dominated by theory. This completes our review of the CKM matrix ele- 
ments Vij. 



4. Radiative and Semileptonic Rare B Decays 

Two inclusive rare B-decays of current experimental interest are B — > X s "f 
and B — > X s l + l~, where X s is any charmless hadronic state with the 
strangeness quantum number s = 1. They probe the SM in the elec- 
troweak b — > s penguin sector. The CKM-suppressed decays B — » Xd"f 
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and B — > Xdl + l~ are difficult to measure due to low rates and formidable 
backgrounds. Instead, the search for B — ► decay has been carried 

out in the exclusive decay modes B — > (p, w)7- Combined with the decay 
£> — > K*j, these decays provide constraints on the CKM parameters. The 
CKM-suppressed decays B — > (7r, p, w)£ + £~ have not yet been measured. 
We review some of these rare -B-decays in the context of the SM. 



4.1. B — > X s ~/: SM vs. Experiments 

The effective Lagrangian for the decays B — > X s ~f obtained by integrating 
out the top quark and the heavy electroweak bosons reads as follows in the 
SM: 

AC 8 

£c ff = -^Ofc^CiMQi • (23) 

* i— 1 

In writing this, unitarity of the CKM matrix has been used and the term 
proportional to the small matrix elements V* s V u b has been ignored. The 
complete list of operators and their Wilson coefficients in the NNLO ap- 
proximation evaluated at the scale fi = rrib can be seen elsewhere . 
The dominant four-quark operators Qi and Q2 are 

Qi = {sc)v- A (cb)v-A, Q2 = (siCj) V -A {c-jb^v-A, (24) 

and the electromagnetic and chromomagnetic penguin operators Q7 and 
Qs are 

em b (p) grrhip) 

<?7 = ( S(T [1 + 75] , Qs = ^—( Sa [1 + 75JT 6)1 



The factor m&(^) is the MS mass of the 6 quark. QCD-improved calcula- 
tions in the effective theory require three steps (for a review see ^): 

(i) Matching Ci(po) (Mo ~ Mw,m t ): They have been calculated up to 
three loops 78 . The three-loop matching is found to have less than 2% 
effect on B(B -> X s7 ) 78 . 

(ii) Operator mixing: This involves calculation of the anomalous dimen- 
sion matrix, which is expanded in a s (fi). The anomalous dimensions up to 
a?(/j) are known ^ since a decade, and the a^(p) calculations have been 

on 01 00 

completed recently in a series of papers ■ ' . 

(iii) Matrix elements {Oi){jib) (p& ~ m^): The first two terms in the ex- 
pansion in are known since long . Exact results to Olotz) were 
obtained for Q7 in ■ 00 and for Q 8 in . For Qi the virtual corrections 
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at G(a s ) were calculated in ■ °'' ° , but those at 0{a 2 s ) are known only 
in the large- (3o limit . A calculation that goes beyond this approximation 
by employing an interpolation in the charm quark mass m c was reported in 
90 ; and has been used in estimating the NNLO branching fraction for the 
inclusive decay B — > X s "/ ^1 Finally, one has to add the Bremsstrahlung 
contribution b — > sjg to get the complete decay rate, which in 0(a s ) was 
done in ^ 2 and in 0{a 2 s ) in ^3. 

In the MS scheme, the NNLO branching ratio for E 1 > 1.6 GeV is 
calculated as 

B(B -> X s -y) SM = (3.15 ± 0.23) x 1(T 4 . (26) 

This amounts to a theoretical precision of about 8%, comparable to the 
current experimental precision ^ 

B(B -» X s7 ) Ex pt. = (3.55 ± 0.24lg;^ ± 0.03) x 10~ 4 . (27) 

A comparison of the two shows that the SM estimate is in agreement with 
data though the SM central value lies below the experiment by about la. 
This allows for speculations about a beyond-the-SM contribution interfering 
constructively with the SM amplitude. A case in point is a 2Higgs doublet 
model (2HDM); the preferred value is m#+ ~ 650 GeV with a 95% C.L. 
lower bound m#+ > 295 GeV. However, more conservatively, the proximity 
of B(B — > X s j) in the SM and experiment puts bounds on the parameters 
characterizing new physics. This has been worked out, together with other 
constraints, in the context of supersymmetry ^4 

The current (NNLO) theoretical precision on B(B — > X s j) has also been 
investigated in the context of SCET using a multi-scale OPE involving three 
low energy scales: rrib, \A^&A and A = mi, — 2Eq, where Eq is the lower 
cut on the photon energy. Large logarithms associated with these scales are 
summed at NLL order. The sensitivity to the scale A introduces additional 
uncertainties. Taking this into account, Becher and Neubert ^ estimate 
B(B —> X s j) = (2.98 ± 0.26) x 10~ 4 , which increases the departure of the 
SM from data to about 1.4a. 



4.2. B — > X s l+l~: SM vs. Experiments 

To study the decays B — > X s £ + £~ , one has to extend the operator basis in 
the effective Lagrangian (23) by adding two semileptonic operators 

2 2 

^ = ^(s^b)v-A (ti"t) , Qw = -^(s^b) V -A (£^15 I) ■ 
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The corresponding Wilson coefficients C$ (/i) and Cio (/i) have the following 
perturbative expansion: 



c w = c<°> + aAMElc-S' + ... . (28) 

After an expansion in a s , the term Cg _1 ^(/x) reproduces the dominant part 
of the electroweak logarithms that originate from photonic penguins with 
charm quark loops: 

_^_ C (-.) (mi )^,»^ + OK), (29) 

leading to ■^t^C { f 1) {m b ) ~ 2. With C^\m b ) ~ 2.2, one needs to 

calculate in the NNLO accuracy. The NNLO calculation of the decay 

B — ► X s l + l~ corresponds to the NLO calculation of B — ► A s 7, as far 

as the number of loops in the diagrams is concerned. 

The process B X s £ + l~ differs greatly from the radiative decay 

B — > X s 7 as far as non-perturbative contributions are concerned. The 

largest effect in B — > A s £ + ^ _ from the intermediate cc states comes from 

the resonances J lib, ib' and tb" decaying to £ + £~, which can be either 

Q7 

modeled, for example, as done by Krugcr and Sehgal y ' using dispersion 
relations and data on <j(e + e~ — > cc — > hadrons), or else experimental cuts 
are imposed on g 2 to remove the resonant regions and the short-distance 
contribution is extrapolated through these cuts. Then, there are factoriz- 
able l/m c and 1/mb power corrections, similar to those in B — > X s -f, which 
can be calculated using the OPE and HQET. As is the case for B — > A„7 
and _B — » X u lvi, there are no l/m& corrections. The 0(l/m 2 ) corrections 
in this framework were calculated first in ^8 an d corrected in The 
0(l/ml) corrections were calculated in 100 The l/m c factorizable power 
corrections were calculated in 

Including the leading power corrections in l/m^ and l/m c and taking 
into account various parametric uncertainties, the branching ratios for the 
decays B -» X s £+r in NNLO arc 102 : 

B(-B -> A s e+e-)sM B(S -» X 5/ u+ M ~)sm - (4.2 ± 0.7) x 1(T 6 , (30) 

where a dilepton invariant mass cut, m« > 0.2 GeV, has been assumed for 
comparison with data given below. These estimates make use of the NNLO 
calculation by Asatryan et al. restricted to s = q 2 /ml < 0.25. The 
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spectrum for s > 0.25 has been obtained from the NLO calculations using 
the scale fib ~ mb/2, as this choice of the scale reduces the NNLO contri- 
butions. Subsequent NNLO calculations covered the entire dilepton mass 
spectrum and are numerically in agreement with this procedure, yielding 
B(B -» X s fi+fi-) SM = (4.6 ± 0.8) x 10~ 6 104 < 105 . The difference in the 
central values in these results and (30) is of parametric: origin. 

The BABAR and BELLE collaborations have measured the invariant 
dilepton and hadron mass spectra in B —* X s £ + £~. Using the SM-based 
calculations to extrapolate through the cut-regions, the current averages of 
the branching ratios are ^ : 

B(B -> X s e + e~) = (4.7 ± 1.3) x 10" 6 , 
B(B -» X sf i+fi-) = (4.3±1;1) x 10~ 6 , 
B(B -» XJ+£-) = (4.5tl;g?) x !0~ 6 - (31) 

Thus, within the current experimental accuracy, which is typically 25%, 
data and the SM agree with each other in the b — > s electroweak pen- 
guins. The low g 2 -region (say, q 2 < 8 GeV 2 ), which allows the most precise 
comparison with the SM, suffers both from the statistics and a cut on the 
invariant hadronic mass recoiling against the dilepton. A cut mx > 2GeV 
and mx > 1 , 8 GeV have been used by the BELLE and BABAR collab- 
orations, respectively. The effects of these cuts have been studied in the 
Fermi-motion model which has been used in the experimental analysis 
of the data so far. Subsequently, the B — > X s £ + £~ rate with an mx cut 
in the low-q 2 region has been calculated using the B — > X s j shape func- 
tion 107 This work, whose impact on the analysis of the B — > X s £ + £~ 
data has yet to be studied, reduces some of the theoretical errors in the 
SM estimates given in eq. (30). In the same vein, it has also been recently 
argued ^at the non-pcrturbative uncertainties in the largc-g 2 region 
(q 2 > 14 GeV 2 ) can be significantly reduced by normalizing the partial 
decay width of B — ► X s £ + £~ with the corresponding partial width of the 
decay B — ► X u £v^. With more data from the B factories, these theoretical 
developments will enable a more precise test of the SM in the B — ► X s £ + £~ 
decays. 

The measurements (27) and (31) provide valuable constraints on 
beyond-the-SM physics scenarios. Following the earlier analysis to deter- 
mine the Wilson coefficients in b — > s transitions ^9, 102, 110^ ^ j^g been 
recently argued ^ that data now disfavor solutions in which the coefficient 
C| ff is similar in magnitude but opposite in sign to the SM coefficient. 

Exclusive decays B — > (K, K*)£ + £~ (P 1 = e ± ,fi ± ) have also been mea- 
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sured by the BABAR and BELLE collaborations, and the current world 

7 

averages of the branching ratios are : 

B(B -» K£+£~) = (3.9 ±0.6) x 10~ 7 , 
B{B^K*e + e~)= (11.3±|:|) x 10" 7 , 
B(B -» K*ft + fT) = (10.3±|;|) x 10~ 7 , (32) 
B(B -» K*£+r) = (9.4t*; 7 ) x 10~ 7 . (33) 

They are also in agreement with the SM-based estimates of the same. A 

119 

calculation based on the light cone QCD sum rules for the form factors llz 
yields 102 : B(B -> K£+£~) = (3.5 ± 1.2) x 10~ 7 , B(B -> K*e+e~) = 
(15.8 ± 4.9) x 10~ 7 , and B(B -» K*fj+fj,-) = (11.9 ± 3.9) x 10~ 7 with the 
errors dominated by uncertainties on the form factors. In the future, these 
errors can be reduced by using the data on B — > (tt, p)£vi to determine 
the B — > (7r, p) form factors. This information can be combined with es- 
timates of the SU(3)-symmetry breaking to determine the B — > (if, if*) 
form factors, enabling to predict the FCNC decay rates and spectra more 
precisely. For the low invariant mass of the dileptons, say q 2 < 8 GeV 2 , the 
SCET framework can be employed to reduce the number of form factors 
and improve the perturbative aspects of these decays. 

The Forward-Backward (FB) asymmetry in the decay B — > X s £ + £r , 
119 

defined as 

j , 2, 1 f 1 , n d 2 B(B ^ X s £+£-) . n . 

AFBiq } = dB{B^X a l+l-)ldq* Jj C ° S9e dcf dco S 6 e Sgn{C ° s6e) 

provides additional constraints on the Wilson coefficients. In particular, 
the location of the zero-point of this asymmetry (called below q 2 ) is a 
precision tests of the SM. In NNLO, one has the following predictions for the 
inclusive decays B -» XJ+£r: gg = (3.90 ±0.25) GeV 2 [(3.76±0.22 th cory± 
0.24 m J GeV 2 ], obtained by Ghinculov et al. 113 [Asatrian et al. 114 ]. 

In the SM (and its extensions in which the operator basis remains 
unchanged), the FB-asymmetry in B — > K£ + £~ is zero and in B — > 
K*£ + £~ it depends on the decay form factors. Model-dependent stud- 
ies yield small form factor-related uncertainties in the zero-point of the 
asymmetry sq = q^jm^ 1^ HQET provides a symmetry argument why 
the uncertainty in so can be expected to be small which is determined 
by 116 Cg // (s ) = -TT^t C 7 ff ■ However, 0(a s ) corrections to the HQET- 

symmetry relations lead to substantial change in the profile of the FB- 

117 118 

asymmetry function as well as a significant shift in sq ■ . They have 
been worked out for B — > K*£ + £~ using SCET 119 . Restricting ourselves 
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to the kinematic region where the light K* meson moves fast and can be 
viewed approximately as a collinear particle, a factorization formula for 
the decay amplitude of B — > K*£ + £~ , to leading power in l/ra&, has been 

derived in SCET ^® . This coincides formally with the formula obtained 

117 1 90 191 

earlier by Bcnekc ct al. , using the QCD factorization approach ■ , 

but is valid to all orders of a s : 

(K*J+£-\H eff \B) = T a V)Ca(<Z 2 ) + (34) 

+ E / — ^±M / du^ K ,{u)T» ± {w,u,q 2 ) , 
± Jo u Jo 

where a =||, _L denotes the polarization of the K* meson. The functions 
and T^{j_ are perturbatively calculable. Ca(<Z 2 ) are the soft form factors de- 
fined in SCET while <I>±(uj) and <j> a K , (u) are the light -cone distribution am- 
plitudes (LCDAs) for the B and K* mesons, respectively. In particular, the 
location of the zero of the forward-backward asymmetry in B — > K*£ + £~ , 
Qq, can be predicted more precisely in SCET due to the improved theoretical 
precision on the scale dependence of qfi. 

Including the order a s corrections, the analysis in estimates the 
zero-point of the FB asymmetry to be 

Qo = (4.07±°;^) GeV 2 , (35) 

of which the scale-related uncertainty is A(g 2 ) sca i c =io!o5 

GeV 2 for the 

range m&/ '2 < Hh < 1m,b together with the jet function scale fii = 
y/JIh x 0.5 GeV. This is to be compared with the result given in also 
obtained in the absence of 1/toj corrections: q^ = (4.39+ ' 35 ) GeV 2 . Of this 
the largest single uncertainty (about ±0.25 GeV 2 ) is attributed to the scale 
dependence. The difference in the estimates of the scale dependence of q 2 
in H9 an( j 118 j s both due to the incorporation of the SCET logarithmic re- 
summation (done in and the different (scheme-dependent) definitions 
of the effective form factors for the SCET currents used in these references. 
Power corrections in 1/rrib are probably comparable to the 0(a s ) correc- 
tions, as argued in So far, g 2 has not been measured experimentally. 
BELLE has published the first measurements ^> ^2 f tne forward- 

119 

backward asymmetry (FBA) . The best- fit results by BELLE for the 
Wilson coefficient ratios for negative value of CV, ^ = — 15.3+4;g± 1-1 an d 

= 10.3t^; 2 ± 1.8, are consistent with the SM values C9/C7 ~ -13.7 and 
C10/C7 — +14.9, evaluated in the NLO approximation. However, for the 
positive value of C7, the measurements lead to ^ = — I6.3I3 7 ± 1.1 and 

= +11.1^3.9 + 1.8 and the two solutions are of comparable significance. 
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With more data at the current B factories, and yet more anticipated at the 
LHC, these measurements are expected to become very precise, providing 
a precision test of the SM in the flavor sector. 

4.3. B — > V~f: SM vs. Experiments 

The decays B — ► V7 (V = K*,p,cu) have been calculated in the NLO 
approximation using the effective Lagrangian given in (23) and its ana- 
logue for b — > d transitions. Two dynamical approaches, namely the QCD 
Factorization ^0 an( j pQQj) 123 have been employed to establish fac- 
torization of the radiative decay amplitudes in the heavy-quark limit. We 
illustrate the QCD-F method, where this factorization is worked out for the 
B — > 124 > 125 < 118 > 126 > 127 > 128 > 129 (see 130 ' 131 for phenomenolog- 
ical updates in NLO, and 10A 100 for the alternative "pcrturbativc QCD" 
approach). In particular, the matrix element of a given operator in the 
effective weak Hamiltonian can be written in the form 

(Vf\Q i \B)=F B ^Tl + J dudu<j>l(u)<fc(u)T}\u,u). (36) 

The non-perturbative effects are contained in F B ^ V± 7 the B — > V transi- 
tion form factor at q 2 = 0, and in 4> B and the leading-twist LCDAs 
of the B- and F-mesons. The hard-scattering kernels T\ and T} 1 include 
only short-distance effects and are calculable in perturbation theory. Con- 
tributions to the kernel T 1 are closely related to the virtual corrections to 
the inclusive decay rate, and are referred to as vertex corrections. Those to 
the kernel T 11 are related to parton exchange with the light quark in the 
_B-meson, a mechanism commonly referred to as hard spectator scattering. 
It is expected that the factorization formula is valid up to corrections of 
0(A QC u/m b ). 

The derivation of the factorization formula from a two-step matching 
procedure in SCET has provided additional insight into its structure. The 
technical details for B — > V"/ in NLO have been provided in ^ 34: 135 j n 
the SCET approach the factorization formula is written as 

(Vj \Qi\ B) = A t C A ( v± + J doodu^^cPlWtY^u), 

(37) 

where F and fy ± are meson decay constants. The SCET form factor (y ± is 
related to the QCD form factor through perturbative and power corrections 
117, 136, 137, 138, 139, 140, 141 i n SCET the perturbative hard-scattering 

kernels are the matching coefficients AiC A and t\ l . They are known com- 
pletely to next-to-leading order (NLO) (0(a s )) in renormalization-group 
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(RG) improved perturbation theory 1 . Recently, important steps to- 
wards a complete analysis at the next-to-next-to- leading order (NNLO) 

1 49 

in B — > V7 decays have been derived in by obtaining full results for 
the hard-scattering kernels for the dipole operators Q7 and Q s , and partial 
results for Q\, valid in the large- /3q limit and neglecting NNLO corrections 
from spectator scattering. In addition, this work provides the virtual cor- 
rections to this order for the B — ► Vj decays, as they can not be obtained 
from the published calculations for the inclusive decay B — > X s -f, discussed 
previously. 

In SCET the hard-scattering kernel t] 1 for a given operator is sub- 
factorized into the convolution of a hard-coefficient function with a uni- 
versal jet function, in the form 

t} I (u,w)= ( dT&iC B1 (T)j±(T,U,Lj) = &iC B1 *j ± . (38) 

Jo 

The hard coefficients AiC B1 contain physics at the hard scale nib, while 
the jet function j± contains physics at the hard-collinear scale \fm~bK. The 
hard coefficient is identified in a first step of matching QCD — > SCETi, and 
the jet function in a second step of matching SCETi — > SCETn- Details in 
NLO have been worked out for B — > V"/ in l^ 4. 

The effective field-theory techniques are crucial for providing a field- 
theoretical definition of the objects in (37), and for resumming large per- 
turbative logarithms of the ratio mi, /Aqcd in t\ l . In the effective-theory 
approach resummation is carried out by solving the renormalization-group 
equations for the matching coefficients AiC B1 . Since these coefficients enter 
the factorization formula in a convolution with the jet function j±, their 
anomalous dimension is a distribution in the variables r and u. The evolu- 
tion equations must be solved before performing the convolution with j±. 
Therefore, resummation is not possible in the original QCD factorization 
formula (36), where the hard-scattering kernels Tj l are obtained only after 
this convolution has been carried out. 

Using the SCET framework, the branching ratios in the NNLO are as 
follows 142 : 

£nnlo(B° -> ^*°7) - (4-3 ± 1.4) x 1(T 5 , 
Bnnlo(S ± -> A" ± 7) ~ (4.6 ± 1.4) x 10" 5 , 
Snnlo(B s ° -» 07) - (4-3 ± 1.4) x 10" 5 . 

It should be noted that, very much like the B — ► X s 7 case, the complete 
NNLO calculations for the virtual corrections to the matrix element of 
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the operator 0\ in B — ► i-T*7 are not yet at hand. In addition, the hard 
spectator corrections from this operator are calculated only in NLO. In the 
NNLO branching ratios quoted above, the errors are increased to take these 
missing pieces into account. 

The above theoretical branching ratios, compared with the current ex- 

f\ 1 1/10 

perimental measurements ' , yield the following results : 



gNNL0(g° ^ K**j) 
£Expt(5 - K*°j) 
gNNL (i3 ± ^ K^j) 

B Expt (B± K*± 7 ) 

BNNL0(Bg -> 07) 

B E xpt(B a ° -» 07) 



= 1.1 ±0.35 ±0.06, 
= 1.1 ±0.35 ±0.07, 
= 0.8 ±0.2 ±0.3 . 



The decays B — ► (p,uj)j involve in addition to the (short-distance) 
penguin amplitude also significant long-distance contributions, in particu- 
lar in the decays B ± — ► p ± 7. In the factorization approximation, typical 
Annihilation-to-Penguin amplitude ratio is estimated as ^ 4 ^: 6a(/5 ± 7) = 
0.30 ±0.07. 0(a s ) corrections to the annihilation amplitude in B — > ,07 are 
not known; also the proof of factorization of this amplitude is still not at 
hand. The annihilation contribution to the decays B° — * p°j and B° — > 0^7 
is expected to be suppressed (relative to the corresponding amplitude in 
i? ± — > /5 ± 7) due to the electric charges (Qd/Qu = — 1/2) and the color 
factors, and the corresponding A/P ratio for these decays is estimated as 
e A (A) - -ca(w7) ^ 0.05. 

The decay amplitudes for £? — > (p, w)7 depend on the CKM matrix 
elements V t * d V t b (from the penguin diagrams) and V* h V u d (from the an- 
nihilation diagrams). Hence, these decays provide potentially very power- 
ful constraints on the CKM parameters, p and fj. Since a large number 
of observables can be measured in these decays, such as the individual 
branching ratios for B ± — ► p ± j and B° — ► (p°,co)j, isospin- and SU(3)- 
violating asymmetries in the decay rates, and direct and time-dependent 
CP asymmetries, they have been studied theoretically in a number of pa- 
pers 124 ' 125 ' 127 ' 128 ' 130 ' 144 . Experimentally, a beginning has been 
made in the measurements of the b — > d"f transition through the measure- 
ments of the branching ratios for B —* (p, oj)j, reported by BABAR and 
BELLE. Current measurements are not very precise, as can be seen from 
the current world averages 7 (in units of 10~ 6 ): B{B ± — > p ±r y) = 0.96±0.23, 
B(B° -> p° 7 ) = 0.77 ± 0.14 and B(B° -» W7) = 0.41 ± 0.15. In addition, 
theoretical estimates suffer from large hadronic uncertainties, dominated 
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by the imprecise knowledge of the form factors. Hence, the resulting con- 
straints on the CKM parameters are not very quantitative. Theoretical 
uncertainties are greatly reduced in the ratios of the branching ratios in- 
volving the decays B — > (p, 07)7 and B — > K*"f. Calling the ratios of the 
branching ratios R ± (pj/K*j) and R (p"f / K**y), for the decays of the B^ 
and B° mesons, respectively, one has 



R ± { P1 /K* 1 ) 



V td 



i?°( P7 /^*7) = \ 



V td 



(39) 

2 (Ml - M 2 ) 3 , „ n 

where C = T({Q)/T^* (0), with Tf(0) and T/^O) being the transition 
form factors evaluated at q 2 = in the decays B — > p7 and £> — > -ftT*7, 
respectively. The functions AR ± (e A , p, fj) and Ai?°(e^, p, 77), appearing on 
the r.h.s. of the above equations encode both the 0(a s ) contribution to 
the penguin amplitudes and annihilation contributions, with the latter es- 
timated so far only in the lowest order. They have a non-trivial depen- 
dence on the CKM parameters p and f\ 125 _ Theoretical uncertainty in 
the evaluation of the ratios R ± (pj/K*j) and R°(pj/K*j) is dominated 
by the error on the quantity £, and to some extent also by the errors 
on the parameters and e^, characterizing the annihilation/penguin ra- 
tios. In the SU(3) limit £ = 1; SU(3)-breaking corrections have been cal- 
culated in several approaches, including the QCD sum rules and Lattice 
QCD. With the current values for the ratios R ± (p 1 /K* 1 ) = 0.032 ±0.008, 
R°((rf/K*-f) = 0.039 ± 0.007 and R°{uj^/K*^) = 0.021 ± 0.007, the cur- 
rent world average of | V*d/V( S | from the ratio of B — > (p, u;)-f and B — > K*-f 
is 145 : 

\V td /V ts \ = 0.194tHii(exp) ± 0.014(th) , (40) 

where the Light-cone QCD sum rules ^1 have been used to estimate the 
hadronic input quantities. This determination is compatible with the one 
from the mass difference ratio AMs s /AMs d given in cq. (22), but less 
precise. 



5. B — > MiM 2 Decays 

Exclusive non-leptonic decays are the hardest nuts to crack in the theory 
of B-dccays. Basically, there are four different theoretical approaches to 
calculate and/or parameterize the hadronic matrix elements in B — > M1M2 
decays: 
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(1) SU(2)/SU(3) symmetries and phenomenological Ansaetze 
146, 147, 148, 149 

(2) Dynamical approaches based on perturbative QCD, such as the QCD 
Factorization ^® and the competing pQCD approach 

(3) Charming Penguins ^® using the renormalization group invariant 
topological approach of Buras and Silvestrini * . 

(4) Soft Collincar Effective Theory (SCET) 43 < 44 < 45 > 46 , which we have 
already discussed in the context of radiative and semileptonic decays. 

These approaches will be discussed on the example of the B — > nn 
and B — > Kir decays for which now there exist enough data to extract 
the underlying dynamical parameters. Prior to this, however, we discuss 
the measurements of the angle (3 (or cf)i) from the experiments at the B- 
factories. 



5.1. Interplay of Mixing and Decays of B°- and 
B° -Mesons to CP Eigenstates 

We start with the discussion of the transition b — > ccs, which is dominated 
by the tree topology. The time-dependent CP asymmetry in the decays 
B° — > / and B° — > /, where / is a CP eigenstate, such as J/ipK s and 
J/^Kl, is defined as: 

S K ' T[B°(t) -» /] + r[B°(i) -» /] ' 1 y 

The time evolution of the two flavor eigenstates B° and B° is determined 
by (2 x 2) Hermitian matrices M and T. The physical states (with defi- 
nite masses and lifetimes) are the linear combinations of B° and B°, with 
\Bd(L, H)} = p\B°)±q\B°), dependent on two complex parameters p and q. 
Defining the decay amplitudes A(f) = (f\H\B°) and A(f) = (f\H\B°) of 
the B°- and B° -mesons into the final state /, the time-dependent CP asym- 
metry is determined by the quantity A f , involving the interplay of mixing 
and decay amplitudes: 

For the B° d - B° d mixing, the ratio q/p involves the phase (3 (or <pi), which 
is one of the angles of the unitarity triangle: 

q _ V t b V td _ -2^ mixing _ -2i/3 / 4 ox 

P ~v tb v; d - c • (43) 



25 



The time dependent CP asymmetry (41) is then expressed as 

A f (t) = C f cos(AM Bd t) + S f sm{AM Bd t), (44) 

where AMs d = (0.507 ± 0.005) ps -1 is the mass difference between the 
heavy and light .B^-meson mass eigenstates and the difference in the decay 
widths ATB d has been neglected. The quantities Cf and Sf, called the di- 
rect and mixing-induced CP asymmetries, respectively arc defined in terms 
of the complex variable A/ as follows: 

_ 1-|A/| 2 21mA, 

C/ -TT^7F' i+1a7F (45) 

If the decays B° — > / and B a — > / are dominated by a single amplitude, the 
ratio p(f ) = r}f e _2 *^ deca y, where r\f = ±1 is the CP parity of the state /, 
is a pure phase factor and the asymmetries (45) reduce to the expressions: 

Cf = 0, Sf = -r]f sin 2(0 mixing + 0docay)- (46) 

The decays B° /B° — ► J /ipK s , J/i/jKl, and a number of related final states 
with / being ip(2S)K s , rj c K s , XdK s , and J/iPK*°(K*° i^Tr ) belong to 
the category of plated decays . In all these modes, the direct CP 
asymmetry C/, to a very high accuracy, vanishes, and the quantity S/, the 
mixing-induced CP asymmetry, measures sin(2/3). Averaging over all the 
decay channels, the results of the BABAR and BELLE measurements are 
as follows 4; 

C = 0.049 ±0.022 ±0.017, S = 0.714 ± 0.032 ± 0.018, [BABAR] (47) 
C= -0.019 ± 0.025, S = 0.651 ± 0.034. [BELLE] (48) 

In the BABAR result, the first error is statistical and the second is system- 
atic, while in the BELLE data both the errors have been combined. The 
current world average for Sf — sin(2/3) for the quark transition b — > ccs 
is 7 : 

sin(2/3) = 0.681 ± 0.025, (49) 

where the data from LEP and Tevatron have also been included. Restricting 
(3 in the range < (3 < n/2, two possible values can be extracted (3 = (21. 5± 
1.0)° and (3 — (68.5 ± 1.0)°. The two-fold ambiguity has now been resolved 
by several cos(2/3) measurements, involving the Dalitz analysis of the decay 
modes B° -» D°_ body h°, B° -» K.tt+tt-, B a -» K S K + K~ , and the older 
results on B° -» J/i/jK*°, leading to the determination (3 = (21.5 ± 1.0)°. 

The direct measurement of sin(2/3) in eq. (49) is to be compared with the 
indirect estimate of the same, obtained from the fits of the CKM unitarity 
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triangle (UT). For this, the UTfit collaboration bt) quotes sin(2/3) = 0.739± 
0.044, obtained from the sides of the UT alone, and sin(2/3) = 0.736 ± 
0.042, by including also the CP-violating quantity in the if-decays. The 
results from the CKMfitter group ^ are similar. Thus, SM passes this test 
comfortably. 

Another key test of the SM in the flavor sector is to compare the 
CP-violating quantities Sf and Cf involving the penguin-topology dom- 
inated quark transitions b — > sss and b — > sdd with the ones from 
the transition b — > ccs, dominated by the tree topology and discussed 
quantitatively above. The point here is that penguin amplitudes may re- 
ceive contributions from New Physics. For example, new phases, present 
generically in supersymmetric theories, may reveal themselves, leading to 
Sf =C cs 7^ Sf=sss;sdd- Examples of the final states induced by the transi- 
tion b — ► sss are (0, 77, r/ , KK)K S , and the ones induced by the transi- 
tion b — > sdd are K °(7r°, p°, u) (B° — > r/'K s receives contributions from 
both the transitions). The current measurements of Sf = — r\f sin(2/3 efI ) 
from the penguin-dominated decays are : S^k = 0.39 ± 0.17, S^k = 
0.61±0.07, S Ks k s k s = 0.58±0.20, S n o Ka = 0.38±0.19, S p o Ks = 0.61±g;|^, 
S uK , = 0.48 ±0.24, S foK o = 0.84 ±0.07, S K+K - K o = 0.73 ± 0.10, and 
S v o^o Kg = —0.52 ±0.41. These measurements are not as precise as the ones 
from the b — > ccs decays due to the much smaller branching ratios (typically 
10~ 5 ) compared to the decay B° — > J/ipK s . Also, they involve more than 
one decay topologies. However, within (large) errors, the values of sin(2/3 efI ) 
from the penguin-dominated transitions are consistent with the value of 
sin(2/?) from the tree-dominated transition given above in eq.(49), with 
the possible exception of 5^o , which deviates by about 2tr, and the poorly 
measured odd-man out 5 , 7r o 7r o^ s . It seems that the fog on Sf in the penguin- 
dominated decays from the initial epoch of the B-factory experiments has 
largely evaporated, and the emerging contours of CP asymmetries in these 
decays are very much the same as predicted by the SM. 

5.2. B — > tztz: SM vs. Experiments 

The determination of the phase a is based on the branching ratios and 
CP asymmetries in the quark transition b — > uud. They metamorphise in 
the decays B — > tttt, B — > pn and B — > pp, apart from other final states. 
We concentrate here on the decay B — > -kit, which has received a lot of 
theoretical attention. 

There are three dominant topologies in the B — > 7T7t decays termed as 
Tree (T), Penguin (P) and Color-suppressed (C). In addition, there are 
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several other subdominant topologies which will be neglected in the discus- 
sion below. Parametrization of the T, P, and C amplitudes is convention- 
dependent. In the Gronau-Rosner c-convention these amplitudes can 
be represented as 

V2A+ = -\T\e iST e n [l + \C/T\e lA ] , 

A+- = -\T\ e l&T [e^ + \P/T\ e tS ] , (50) 
V2A 00 = -\T\ c t5T [\C/T\ e lA e* 7 - \P/T\ e* 5 ] . 

The charged-conjugate amplitudes A y differ by the replacement 7 — > —7. 
The amplitudes (50) and the charged-conjugate ones obey the isospin rela- 
tions: 

A+° = -^A+- +A m , A~° = -^A+- +A m . (51) 

There are 5 dynamical parameters \T\, r = \P/T\, 5, \C/T\, A, with 
St = assumed for the overall phase. Thus, the weak phase 7 can be 
extracted together with other quantities if the complete set of experimental 
data on B — > tttt decays is available. 



Tabic 1. Branching ratios (in units of 10 6 ) and CP asymmetries in the 



B — * rnr decays 


B(B+ -V 7T+7T ) = 


- 0.0»_Q 4Q 




= 


= 0.06 ± 0.05 


B(B° -> 7r°7r°) = 


1.31 ±0.21 


A CP (n 


7T°) = 


0.48l°;g 


B(B° -> 7T+7T-) 


= 5.16 ± 0.22 


C C p(^ 




= -0.38 ±0.07 






Scp(tt + 


T") = 


= -0.61 ±0.08 



The experimental branching ratios and the direct CP asymmetries 
Acp(tt°tt°) and Ccp(7r + 7r - ), as well as the value of the coefficient 
S'cp( 7r+7r ~) in time-dependent CP asymmetry, presented in Table 1, have 
been fitted to determine the various parameters (the direct CP asymme- 
try Acp(n + Tr a ) is not relevant for this analysis but can be important in 
determining the size of electroweak contribution in the decays considered) . 
An updated analysis by Parkhomenko based on the paper yields the 
following values for the hadronic parameters: 

\P/T\ = 0.473t°:°°°, 6 = (-40.2^)° , (52) 
\C/T\ = 0.966t°°6i, A = (-56.3t|:*)°, 
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and for the CKM unitarity triangle angle 7 (or equivalently a) 

7 = (65.91*;°)° , a = n - 7 = (92.6^)° . (53) 

Similar fits based on their data have been performed by the BABAR and 
BELLE collaborations resulting in slightly larger values: a = (96tg°)° 
(BABAR) and a = (97 ± 11)° (BELLE). The overall fits performed by 
the CKM-Fitter and UT-Fit groups prefer slightly smaller values, yielding: 
a = (90.7^;|)° 67 and a = (88.7 ± 6.2)° 66 , respectively. All the above 
estimates are in good agreement with each other within the quoted errors, 
stating that the data on B — > nn (as well as the other decay modes B — ► pn 
and B — > pp) arc in agreement with the indirect estimate of the phase a 
from the unitarity triangle. The strong phases 8 and A in Eq.(52) come out 
rather large. In particular, they are much larger than the predictions of the 
QCD-F approach with pQCD in better agreement with data, but 
neither of these approaches provides a good fit of the entire B — > 7T7t data. 

Data on B — > nn decays are in agreement with the phenomenological 
approach of the so-called charming penguins an d with the SCET-based 
analyses by Bauer et al. which also attributes a dominant role 

to the charming penguin amplitude. However, a proof of the factorization 
of the charming penguin amplitude in the SCET approach remains to be 
provided. In addition, SCET makes a number of predictions in the B — > irir 
sector, such as the branching ratio B(B° — ► 7r°7r°): 155 

B(B° -> tt°tt°) = (1.3 ±0.6) x TO" 6 . (54) 

7=64° 

In contrast, predictions of the QCD-F and pQCD approaches are rather 
similar: B(B° — > 7r°7r°) <~ 0.3 x 10~ 6 , in substantial disagreement with the 
data. 

5.3. Present bounds on the phase 7 from B decays 

The classic method for determining the phase 7 (or ^3) 158, 159, 160 
involves the interference of the tree amplitudes b — > uW" — > ucs leading 
to B~ -» D°K- and b -> cW- -» cus leading to B~ -» WE'. These 
amplitudes can interfere if D° and D° decay into a common hadronic final 
state. Noting that the CP= ±1 eigenstates D± are linear combinations of 
the D° and W> states: D° ± = (D° ± ~W)/\/2, both branches lead to the 
same final states B~ — > D±K~ . So, the condition of CP interferometry is 
fulfilled. The decays B~ — > D^K^ are described by the amplitudes: 

A{B~ -» D^X") = 4= -» i' ^ - ) ± A(B- -» ^Jf-)1 . (55) 

V2 L J 
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A+ = 



Since, the weak phase of the b — > u transition is 7 but the b — > c transition 
has no phase, a measurement of the CP asymmetry through the interference 
of these two amplitudes yields 7. The four equations that will be used to 
extract 7 are: 

B(B~ -» D°±K-) + B(B+ -» D%K+) _ 
^ S BF^F)TIFTD^) = 1 + -dk ± 2r DK cos S DK cos 7 , 

(56) 

B(B- -> £>^~) - ff(g+ -> D|g+) _ ±2r DK sin<5 DK sin 7 

Z^if-) +fi(B+ D%K+) ~ l + rl K ±2r BK cos^K cos 7 ' 

Here, tdk is the ratio of the two tree amplitudes tdk = \T\lT2\ ~ 
(0.1 - 0.2), with Ti and T 2 being the CKM suppressed (b -» u) and CKM 
allowed (6 — ► c) amplitudes, respectively, and 5uk is the relative strong 
phase between them. The construction of the final states involves flavor 
and CP-tagging of the various D° states, which can be done, for example, 
through the decays D° + -> tt+tt", D°_ -» i^ S 7r°, and D° -> K~tt+. With 
three unknowns (r DK , <5dk, 7)7 but four quantities which will be measured, 
i?± and A±, one has, in principle, an over constrained system. 

Experimentally, the quantities R± are measured through the ratios: 

With all three quantities R(K/n) and R(K/tt)± measured, one can 
determine R± = R(K/n)±/ R(K/n) . More useful decay modes to construct 
the B — ► ZJ/T triangle can be added to reduce the statistical errors. Along 
these lines, Atwood and Soni have advocated to also include the decays 
of the vector states in the analysis, such as B~ — ► K*~D°, B~ — ► K~D*°, 
and S- -» K*-D*°, making use of the L>*° -> L>°7 and £>*° -» D°7r° 
modes. 

Present measurements in the B — > Dif and Z? — > Z?7r decays by the 
BABAR and BELLE collaborations yielding R± and A± for the DcpK~ 
mode are summarized by HFAG ^ : 

R+ = 1.09 ± 0.09 , A+ = 0.26 ±0.08 [BELLE, BABAR] , 
R- = 0.90 ±0.10, A_ = -0.16 ±0.09 [BELLE]. 

The corresponding quantities for the Dq V K~ and DcpK*~ are also given 
by HFAG 1 . 

A modification of the Gronau-London-Wyler (GLW) method described 
above has been suggested by Atwood, Dunietz and Soni (ADS), where 
B~ — > D°K~ with Z) — > if + 7r~ (or similar) and the charge conjugate 
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decays are implied. BABAR and BELLE use the following definitions for 
the quantities called Aads and .Rads, (the decay modes B~ — ► D°K~ 
followed by D° — » K + tt~ are used to exemplify the method) 



= B(B- -> [Jf+7r-] c Jf-) + £(£+ -> 
ADS - -» [ir-TT+jc^-) +B(B+ -» [if+7r-] D ^+) 

_ -» [Jf+7r-] D A-) - -» [ J Ff- 7 r+]z 3 ir+) 

^ADS = 



(58) 



-» [A'+7r-] r ,A'-) + fi(B+ -» [A--7r+] D A'+) ' 

The current measurements of these observables are summarized by HFAG . 
In the analysis of data, usually the GLW and ADS methods are combined 
and a x 2 -fit is done to determine the profile of the phase 7. 

A variant of the B — > DK method of measuring 7 is to use the de- 
cays B ± — ► DK^ followed by multi-body decays of the D-meson, such 
as D° -> K S TT-7r+, D° -> K S K-R+ and D° -» if S 7r-7r+7r , in which 
a binned Dalitz plot analysis of the decays D°/D° — > AT,g7r~7r + was pro- 
posed 163 Assuming no CP asymmetry in £>° decays, the amplitude 
of the B+ -» L>°.fi:+ -» (XsTr+Ti--)/^ can be written as 

M+ = f(m%,m 2 _) + r DK e^ +s ^f(m 2 _,m 2 + ) , (59) 

where m 2 ^ and m_ are the squared invariant masses of the Kgir + and Kgit~ 
combinations in the D° decay, and / is the complex amplitude of the decay 
D Q — > /f57r + 7r _ . The quantities r.Djc and Sdk are the relative magnitudes 
and strong phases of the two amplitudes, already discussed earlier. The 
amplitude for the charge conjugate B~ decay is 

M_ = f(m 2 _,m 2 + ) + r DK c l( -' 1+5DK) f(m 2 +7 m 2 _) . (60) 

Once the functional form of / is fixed by a choice of a model for D a — ► 
Ksir + n~ decay, the Dalitz distribution for £?+ and B~ decays can be fit- 
ted simultaneously by the expressions for M + and M_, with tdk, 5dk and 
7 (or (^3) as free parameters. The model-dependence could be removed by 
a binned Dalitz distribution _ This is usually called the GGSZ method, 
and has been used to determine 7. The combined fit of both of these meth- 
ods by CKMfitter 67 yields 7 = (76.8tlti)°, to b e compared with their 
overall fit from the CKM unitarity 7 = (67.6*4 5)°- The corresponding fit 
by the UTfit group 66 yields 7 = (67 ± 7)°. Thus, we see that within the 
current experimental error of the direct measurements, which is quite large, 
also the phase 7 is compatible with its indirect estimates in the SM. The 
experimental precision will greatly improve at the LHC, in particular, by 
using 2-body _B s -decays. 
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5.4. B — > Kir: SM vs. Experiments 

We now discuss the decays B — > _K"7r. First, we note that the direct CP- 
asymmetry in the £> — > if 7r decays has now been measured by the BABAR, 
BELLE and CDF collaborations: 

r (— 10.7± 1.8±g;X)% [BABAR], 
A C p(n + K-) = I (-9.3 ± 1.8 ±0.8)% [BELLE], (61) 
[ (-8.6 ±2.3 ±0.9)% [CDF], 

to be compared with the predictions of the two factorization-based ap- 
proaches: A C p(7r+ K~) = (-12.9-=-- 21.9)%[pQCD] 123 and A C p(tt+K-) = 
(-5.4 -=- +13.6)%[QCD - F] 120 , with the latter falling short of a satisfac- 
tory description of data. 

The charged and neutral B — > irK decays have received a lot of theoret- 
ical attention. In particular, many ratios involving these decays have been 
proposed to test the SM 165 < 166 < 167 < 168 and extract useful bounds on 
the angle 7, starting from the Fleischer-Mannel bound -^5. 

The current experimental average R = 0.899 ± 0.049 allows to put a bound: 
7 < 92° (at 95% C.L.). This is in agreement with the determination of 7 
from the B — > 7T7T and B -» decays given earlier and the indirect 

unitarity constraints. Thus, both R and Acp^+if - ) are in agreement with 
the SM. The same is the situation with the Lipkin sum rule 167. 

implying significant electroweak penguin contribution in case i?L deviates 
significantly from unit. With the current experimental average Rl = 1.071± 
0.049, this is obviously not the case. This leaves then the two other ratios R c 
and R n involving the B — ► irK decays of and B° mesons: 



B(B ± -> 7r°if ±) _ 1 ff(ff° -> tt^f K ± ) 

B(B ± ^ 7r±^°) ' " = 2 £(5° 7r°iC ) 



Their experimental values i? c = 1.117 ± 0.071 and i?„ = 0.980 ± 0.067 are 
to be compared with the current SM-based estimates ^4 r c _ 1.14 ± 0.05 
and i?„ = l.ll±g;g|. This implies fl c (SM) - i? c (Exp) = 0.02 ± 0.09 and 
R n (SM) - i?„(Exp) = 0.13 ± 0.08. Possible deviations from the SM, if con- 
firmed, would imply new physics, advocated in this context, in particular, 
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by Yoshikawa Beneke and Neubert ^® and Buras et al. How- 
ever, as of now, one has to conclude that SM is in agreement with the 
measurements of both R c and R n . 

Finally, a bound on B(B° K°~K®) based on SU(3) and B —> mr 
data, obtained recently by Fleischer and Rccksiegel yielding B(B a — > 
K°K°) < 1.5 x 10~ 6 is well satisfied by the current world average 7 B(B° — » 
= (0.96lg;S) x 10- 6 . 

6. B° Physics: Eldorado for the Tevatron and the LHC 

The main goal of b physics at the hadron colliders Tevatron and the LHC is 
to chart out the physics of the B® and Bf mesons and of the 6-baryons. The 
current information on the spectroscopic and decay characteristics of these 
hadrons is still very much in the offing, though clearly the two Tevatron 
experiments have made some incisive inroads in these otherwise uncharted 
territories. Despite the overwhelming performance of the B factory exper- 
iments, there still remain a few landmark measurements to be carried out 
involving B^ and B^ mesons. These include, precise measurements of the 
CP asymmetries in the penguin-dominated exclusive decays, quantitative 
determinations of the Wilson coefficients in the effective theory for weak 
decays (CV, C&, Cg, Cio), which will be made possible by the precise mea- 
surements of the radiative and semileptonic decays B — > (X s ,K*)-f and 
B — ► (X S ,K, K*)l + l~ . It is challenging to measure the inclusive decays at 
the LHC, but certainly exclusive decays will be well measured. 

In this section, a brief list of some selected b physics topics to be studied 
at the LHC is given and discussed. 

• B° -~B° Mixing. 

Apart from the precise measurement of AM Bs = (17.77 ± 0.10 ± 0.07) 
(ps) _1 by the CDF collaboration, there are two other quantities still to 
be measured in this complex: Lifetime difference Ar^, and the phase <j> s . 
These quantities have been calculated to a high precision in the SM ^ . A 
recent update of this work yields ^ 

Ar 

A] J S = (49.7 ± 9.4) x 10~ 4 ; Ar Bs = (0.096 ± 0.039)(ps)- 1 , 

4> s = (4.2 ± 1.4) x 10~ 3 = 0.24° ± 0.08° . (65) 

The current measurements of these quantities from the DO collaboration 
are 173 

Ar Bs = (0.12±g;5g ± 0.02)(ps)- 1 (assuming <f>. = 0) , (66) 
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and 

Ar Bs = (0.17 ± 0.09 ± 0.02)(ps) _1 ; <p s = -0.79 ± 0.56±g;^ . (67) 

The corresponding measurement (assuming <j) s = 0) for Ar^ from CDF 
is 174 Ar Bs = (0.076+°;^ ±0.006) (ps)" 1 , where the hrst error is statisti- 

1 9 

cal and the second systematic. At the LHCb , one anticipates a statistical 
sensitivity of cr(sin0 s ) ~ 0.031 and cr(Ar s /r s ) ~ 0.011, assuming an in- 
tegrated luminosity of 2 (fb)" 1 and using the decay B s — > J/tpcj). This 
sensitivity will be improved by accumulating more data and adding the CP 
modes B s — » J/ip(j> an d the pure CP modes B s — > J/ipri and B s — > J/i^rj'. 
The ATLAS and CMS sensitivities on <j> a are expected to be somewhat 
worse by typically a factor 2. Thus, experiments at the LHC will be able 
to test the SM estimates for both the quantities ATb s and <f> a . 

• Precise measurement of the phase 7. 

At the LHC one can combine the methods of determining the phase 7 
already in use at the B factories involving the decays B + — ► D^K + and 
B° -> DK { *\ with the decay B° s -> D S K. In addition, one can use the 

1 

U-spin symmetry arguments advocated, in particular, by Fleischer , to 
combine data from the decays B d — > ir + n~ and £?" — > K + K~ to constrain 
7. An educated guess is that a sensitivity a(j) ~ 4° with 2 (fb) _1 of 
data is reachable at the LHCb, improving to <j(j) ~ 2.4° with 10 (fb) _1 . 
This will result in an order of magnitude improvement over the current 
precision on this phase. Modest improvements are also anticipated for the 
other two phases a and (3 at the LHC. 

• Leptonic decay B® — > fi + n~ 

1 77 

New and improved upper limits have been presented by the CDF and 
DO collaborations for the decays B® — > and £>° — » They 

are as follows (at 95% C.L.) 

-» m + aO < 9.3 [5.8] x 10~ 8 D0[CDF] , 
B(B" d -» < 1.8 x 10~ 8 [CDF]. (68) 

The CDF and DO upper limits have been combined to yield ^ B{B® — ► 
/i + /i~) < 4.5 x 10~ 8 , to be compared with the SM predictions 1™ 
£(5° -» = 3.4 x 10- 9 and B{B° d -> /x+M") = 1-0 x lO" 10 with 

±15% theoretical uncertainty. Hence, currently there is no sensitivity for 
the SM decay rate. However, as the leptonic branching ratios probe the 
Higgs sector in beyond-thc-SM scenarios, such as supersymmetry, and they 
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depend sensitively on tan/3, the Tevatron upper limit on B{B® — > /i + /i~) 
probes the large tan/3 parameter space, though the precise constraints are 
model dependent l 8 ^ At the LHC, the two main collider experiments 
ATLAS and CMS will reach the SM sensitivity, certainly with the higher 
LHC luminosity, Llhc = 10 34 cm~ 2 s _1 , as the decay B® — > remains 
triggerable with the high luminosity. 

• Charmless non-leptonic B s — > h\hi decays. 

The experimental program to study non-leptonic decays B° s — ► hih 2 
has started (here stand for charmless light vector or pseudoscalar 
mesons) with first measurements for the branching ratios B® — > K + ir~~ and 
B® — > K + K~ made available recently by the CDF collaboration ^> 10. Re- 
markably, the first direct CP asymmetry involving the decay B® — * A + 7r~ 
and its CP conjugate mode reported by CDF is found to be large, with 
A CP (B° -> A+7T-) = (39 ± 15 ± 8)%. This large CP asymmetry was 
predicted by Lipkin based on SU(3) symmetry arguments. This al- 
ready tests various dynamical models, such as QCDF 183 , SCET 184 and 
pQCD 185 . With the ongoing fe-physics program at the Tevatron, but, in 
particular, with the onset of the LHC, we expect a wealth of data involving 
the decays of the hitherto less studied B® meson. The charmless B® — ► h\hi 
decays are also important for the CP asymmetry studies and the determi- 
nation of the inner angles of the unitarity triangle. As already stated, a 
number of charmless decays B® — > h\hi can be related to the B^ — > h\h<i 
decays using SU(3) (or U-spin) symmetry, and hence data on these de- 
cays can be combined to test the underlying standard model and search for 
physics beyond the SM under less (dynamical) model-dependent conditions. 
Anticipating the experimental developments, many studies have been de- 
voted to the interesting charmless B® — > h\h,2 decays, waiting to be tested 
at the LHC. 

7. Summary and Outlook 

Summarizing, dedicated experiments carried out over several decades com- 
bined with progress in theoretical techniques embedded in QCD have en- 
abled a precise determination of the CKM matrix elements. The knowledge 
of the third row and the third column of Vckm has come from &-physics, 
which we discussed at length in this review. Of these, precise determina- 
tion of V c b and V u b required good control over the perturbative and non- 
perturbative aspects of QCD. The current precision on the direct deter- 
mination of V t b from the decay t — > bW is limited by statistics and this 
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will vastly improve from the top quark studies at the LHC and later at the 
ILC. The determination of Vt s and Vtd require not only precise knowledge of 
QCD in b decays but implicitly also the assumption of the CKM unitarity, 
as they are determined from the loop-induced b s and b — > d transitions. 
Their current best measurement is through the mass differences AMb s and 
AMe d , and the precision on these matrix elements (typically 10%) is com- 
pletely dominated by theory. A complementary determination of V t d and 
Vts is also at hand from the radiative penguin transitions b — > (d, s)"f and 
the exclusive decays B — > (K* , p,uj)j, but the current precision is limited 
by both experimental statistics and non-perturbative aspects of QCD, such 
as the transition form factors in exclusive decays. This surely will improve 
over the next several years. 

Experiments have also firmly established the phenomenon of CP vi- 
olation in the K and B meson sectors. The various CP asymmetries in 
these decays are found compatible with each other and, with some help 
from QCD, have a consistent interpretation in terms of the single com- 
plex phase of the CKM matrix. Again, in principle, there is ample room 
also for beyond-the-SM weak phases, which would lead to very different 
patterns of CP asymmetries in the tree-dominated versus loop-dominated 
transitions. This has not been borne out by experiments at the B factories. 
While the current data is not equivocal on all the decay channels, and the 
dynamical aspects of not all the measured _B-meson decays are quantita- 
tively understood, the experimental case for the extra weak-phases is rather 
weak. Ongoing experiments at the B factories are expected to significantly 
reduce the errors on the quantities Sf and C/ in penguin-dominated decays 
to settle the issue of new weak phases in B decays. 

From the foregoing one has to tentatively conclude that the CKM 
paradigm is now firmly established as the dominant mechanism of flavor 
transitions in the quark sector. Whether future experiments, such as at the 
LHC and (Super) B factories, will force us to modify this paradigm remains 
to be seen. We expect on theoretical grounds that there is New Physics, 
probably just around the corner, to solve the outstanding issue of the gauge 
hierarchy. The resolution of this problem together with the unification of 
the gauge couplings and the search of viable candidate(s) for dark matter 
requires a TeV scale New Physics. Assuming that supersymmetry is the 
most viable candidate for the impending New Physics to be discovered by 
experiments at the LHC, the central issue in the LHC era would be to pin 
down the underlying flavor aspects of this theory. However, if the current 
experimental trend is any indicator, then very likely the New Physics will 
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be of the minimal flavor violating type, or something akin to it. Verifying 
this and, more importantly, quantifying the subdominant flavor structures 
in supersymmetry (or its alternative) is the next task to which theory and 
experiment have to gear themselves up. 
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